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Abstract 

This is the first of a series of papers devoted to a systematic study of QCD correlation 
functions in a framework of 'instanton vacuum' models. The topic of this paper is to work 
out approximate formulae for quark propagators in a multi-instanton environment. As an 
application, and also as a necessary step toward understanding the correlation functions, 
we study the propagators of scalar and spinor quarks, using the simplest possible model, 
the so called 'random instanton vacuum' (RIV). Results related to heavy-light mesons, 
are found to be very consistent with phenomenology. 
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1. Introduction 



The understanding of non-perturbative phenomena in the QCD vacuum remains to 
be a challenging problem of contemporary physics. One particular way to obtain insights 
in the complicated quark-quark and quark- ant i- quark interactions is provided by point- 
to-point correlation functions, which, at medium distances of 1/3-1 fm, display strong 
deviations from free quark propagation. These deviations are strongly channel dependent 
and change widely in both sign and magnitude (see PJ for a recent review). 

Although the hadronic spectrum can be described more or less in terms of a universal 
confining force and perturbative spin forces, such a simplified picture definitely fails to 
describe the point-to-point correlation functions. The actual situation is much more 
complex, and additional theoretical inputs are definitely needed. 

Tunnelling phenomena in the QCD vacuum, described semiclassically by the so called 
instanton solutions were suggested as explanations of at least some of these effects (or 
even most of them). Classical papers of the one-instanton era describe the discovery 
of the instanton solution its physical interpretation as tunneling and the relation to 
chiral anomalies ||, and the semiclassical integration of the quantum fluctuations 
The first applications of instantons to QCD problems || , based on the so called dilute gas 
approximation, attracted a lot of attention in the late seventies. However, the absence of 
an explanation for the diluteness and the semiclassical nature of the instanton ensemble 
led to such a pessimism that most people left the field around 1980. 

However, at the same time, phenomenological studies of instanton-related effects || [?] 
have shown that the induced effective interaction of light quarks has exactly the proper- 
ties needed to explain several puzzles of the hadronic world. Eventually, this development 
has resulted in the so called 'instanton liquid' model || |9|, and its phenomenological suc- 
cess revived the hopes of describing the QCD vacuum by a dilute system of semiclassical 
fluctuations. 



Multiple attempts [[[(], [Ll|] |L2], [L3| [15] were made to create a quantitative theory 
of the 'instanton component' of the QCD vacuum. The problem appears to be very 
difficult: even if one disregards all other non-perturbative phenomena and considers only 
instantons, it is complicated both by the gauge-induced and, especially, by the fermion- 



induced interaction between them. The theory of these phenomena, which we refer to 
as the interacting instanton approximation (IIA)Q, has been developed during the past 
decade, and the present series paper can be considered as an attempt to clarify its relation 
to experimental information on correlators in a quantitative way. 

So far, the only numerical evaluation of correlators in a multi-instanton environment 



was made in a series of works by one of us [16, 17]. The first results on a set of correlation 



functions have well reproduced the data, sometimes in surprising detail. However, the 
accuracy of the calculations was not high, and both the statistical and the systematic 
errors should be significantly reduced. This is one of the tasks addressed by the present 
series of papers. 



Recently, the first results for point-to-point correlators from lattice simulations [ |18 
were reported. Remarkably, inside uncertainties they generally agree both with experi- 
ment and IIA predictions. Whether these results are or are not related to instantons can 
be clarified by further lattice studies. In the long run much better lattice data will be 
obtained, so that a more detailed comparison with experiment becomes possible. In any 
case, this recent development makes the situation in the field very exciting. 

The present paper is the first in a series of three works. The objective is to perform a 
systematic calculation of all major correlation functions at distances of the order of 1 fm. 
A systematic comparison of the results with experimental data is made whenever possible. 
We also compare our results with predictions of the operator product expansion (OPE), on 
which the so called QCD sum rules are based. As this task happens to be quite formidable 
by itself, merely because of the literally dozens of different correlation functions etc., we 
feel that pushing forward the theory of interacting instantons at the same time would 
be very confusing. Instead, we have adopted the idea to establish a kind of 'benchmark 
calculation' for all correlators, using an as simple input model as possible. The particular 
framework used below will be called the 'random instanton vacuum' (RIV). It is nothing 
else but the 'instanton liquid model', originally suggested by one of us f|, supplemented 
by the simplest assumption of a random distribution of the collective coordinates. We 
do not fit any parameters, leaving them as they were proposed a decade ago at the one- 
instanton level. It is relatively easy to perform high-statistics numerical studies of this 



^or a recent review and references on this subject we refer to jl| 



simple model, and, as we demonstrate below and in the next papers of this series [19, [2C 
it reproduces most of the experimental observations amazingly well. 

The present paper is organized as follows. In section 2 we discuss the main ingredients 
of our model, and we pay particular attention to the zero mode sector of the Dirac 
operator. Numerical results for its spectrum will be presented. The general structure 
of the quark propagator and its numerical evaluation will be discussed in section 3. In 
section 4 we discuss the propagation of a scalar quark in the RIV. The more complicated 
spinor propagator is discussed in section 5, and the results are applied to heavy-light 
systems in section 6. Concluding remarks are made in section 7. 

The mesonic [rj| and baryonic |2(J correlation functions are discussed in two subse- 
quent papers of this series. After such benchmark for all correlators is established, we 
are planning to return to their discussion in the presence of instanton interactions with 
correlations among the positions and orientations. Those studyies are subject of future 
investigations. 



2. The model and the zero mode zone 

The general reason why instantons are so important for the physics of light quarks is 
related to the so called fermionic zero mode, which exists for any topologically non-trivial 
gauge field. In order to explain this let us consider the Euclidean quark propagator defined 
by 

3=-< — 1 . >, (2.1) 
i>D^ + %m 

where iD^ = id^ + is the covariant derivative containing the external gauge field, 
which should be averaged over with the proper weight following from the QCD partition 
function. The operator may be inverted by first diagonalizing it, 

iD^<f> x (x) = \(f)x(x), (2.2) 

which leads to the general expression 

(23) 



for the propagator in the gauge field A^. It is clear that in the chiral limit (m — > 0) the 
eigenvalue spectrum at small A has important effects on the propagation of hadrons. In 



particular, one can easily derive the following general formula |H| for the Euclidean quark 
condensate 

— 7T? 

< W> >e= yp(X = 0), (2.4) 

where p(A) is the spectral density^. 

'Zero modes', discovered in |3[], are solutions of the Dirac equation in the field of 
one instanton with eigenvalue A = 0. Their existence explains the so called U(l) chiral 
anomaly: while tunneling, quarks with one chirality 'dive into the Dirac sea', and those 
with the opposite chirality 'emerge' from it. In [|19|, ^(J we will show that these phenomena 
are responsible for a much wider range of physical effects than was believed previously. 
In particularly, not only the rj' channel is affected, but, say correlators in the ir meson or 
the nucleon channels can be reproduced as well. 

The evaluation of a quark propagator in the multi-instanton background field is a 
complicated problem. In particular, even a more or less dilute set of instantons (and 
anti-instantons) is described by 12 collective variables per instanton (position, size and 
orientation), and the spectrum of the Dirac operator depends on all of them. Roughly 
speaking, the problem resembles the propagation of a sound wave through a turbulent 
atmosphere, containing many vortices and anti-vortices. 

To explain our approximations, one can use an analogy from solid state physics in 
which the role of atoms is played by instantons, and the electron bound state is represented 
by the quark 'zero mode'. The general lesson from condensed matter physics is that at 
finite density of the atoms, atomic bound states may become collective and form 'zones' 
of delocalized states^. Analogously, fermionic states in the 'instanton vacuum' form the 
so called 'zero mode zone' (ZMZ)Q. 

Furthermore, if the ensemble of instantons is dilute, the off-diagonal matrix elements 



2 We do not repeat here the standard arguments on the order of the limits m — ► and V — > oo. For a 
very detailed recent discussion of this topic see |2^] . 

3 It is important that such zones exist as well for disordered systems (liquids), although in this case 
they do not have sharp boundaries. 

4 What is important here, is that the mode for an individual 'well' is exactly at zero, due to some 
topological theorem: therefore all small distortions are irrelevant. We thus deal with a degenerate 
situation, in which the spectrum is governed by the off-diagonal matrix elements. 
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are small and the zone should be narrow. If so, this set of states may dominate the 
spectrum near zero virtuality, and therefore be mainly responsible for p(A = 0) and the 
quark condensate. Direct evidence that this indeed happens in the QCD vacuum has been 



obtained in lattice studies |23|, Much more work in this direction can and should be 
done in order to clarify the relation between lattice results and the theory of instantons. 

Returning to the propagator, one arrives at a comprehensive physical picture of how 
light quarks can easily propagate over large distances in the QCD vacuum: they can do 
it by simply jumping from one instanton to the next. It is analogous to what conduction 
electrons do, in solid or liquid metals. Due to the degeneracy of atomic levels, quite small 
perturbations can crucially enhance or reduce probability of this to happen. 

Continuing this analogy further, we have learned from condensed matter physics that 
going beyond the one-electron approximation one comes across an important interaction 
between conduction electrons due to the possibility, that two of them may happen to 
be at the same atom. The same is true for quarks in the instanton vacuum: the fact 
that they may meet each other at the same instanton leads to their strong interaction. 
Technically this interaction is revealed, when one compares the average propagator squared 
(for mesons) or cubed (for baryons) to the same power of the average propagator. The latter 
is studied in detail below, while the former quantities will discussed in two subsequent 
papers. 

Having outlined the main physics involved, let us now specify the ZMZ in more detail. 
In a basis of zero modes the Dirac operator reduces to the following matrix 



where the overlap matrix elements are defined by 

T H -= [ d^x&ix)^ ■ d(j)j(x). (2.6) 



The ordinary derivative appears because the gauge field in the covariant derivative has 
been eliminated by using the equations of motion for the zero modes. Although it is not 
very important for what follows, in this paper we will use the overlap matrix elements 
corresponding to 'streamline' gauge field configurations. They are discussed in |25| and 



are equal to 

T ^(^U-kA F{X) (2 7) 

TlI ~2 Tl { r Ul Ul ) ( PlPI y^ (2 - 7) 

where R — Ri — Rj is four-distance between the centers, Uj 1 Uj is the relative orientation 
of the instantons and A is a conformally invariant parameter defined by 



1 R 2 + P j + P} l f (R 2 + P j + pW 

2 p lPl 2\ p\p\ 



Here and below we use the standard color matrices = (a, =B) /t . These matrices project 
the relative orientation matrix onto its upper 2x2 block. The scalar function F(X) can 
be reduced to a simple one dimensional integral 

/>oo cLtt'^^ 

FW=6 l (r + l/A)V(r + A)'/» - (2J) 
Asymptotically, for large R the overlap matrix elements behave as T ~ p 2 /R 3 , and thus 
correspond to the exchange of a massless quark. 

In this set of papers we are going to use a simplified instanton ensemble, the so called 
random instanton vacuum (RIV). We assume the following: (i) all instantons have the 
same size po; (ii) they have random positions and orientations; (iii) the instanton and 
anti- instanton densities are both equal to N/2V, where V is the volume of the Euclidean 
space time and iV is the total number of pseudoparticles. 

Thus, there is essentially only one dimensionless parameter in the model / = 7r 2 V/NpQ 
which describes 'diluteness' of the vacuum. All propagators and correlators to be con- 
sidered, normalized with respect to their values for free massless quarks, are also given 
by dimensionless ratios. Those can only depend on /, while the overall distance scale is 
given by a second parameter, e.g. by the 'average separation' R defined by (V/iV) 1 / 4 . 

As suggested a decade ago in ||, the parameters are chosen such that several bulk 
properties of the QCD vacuum are reproduced. In order to obtain the traditional value 
for gluonic condensate one has to take R = 1.0 fm (inside uncertainties), whereas the 
quark condensate (and several other phenomenological parameters) are reproduced for 
po = 1/3 fm. 

Somewhat surprisingly, already this simple model leads to a very reasonable descrip- 
tion of many correlation functions. This does not mean that we have a full description of 
the QCD vacuum, and a number of disclaimers have to be made: 



(i) First of all, the model does not include neither the (so far mysterious) field fluctu- 
ations leading to confinement, nor even perturbative fields, leading to phenomena like the 
Coulomb interaction or radiative corrections to 'asymptotically free' propagation. The 
confining fields are supposed to modify our results substantially for quarks travelling far 
apart, while perturbative corrections are of 0(a s (x)/'K) ~ 10 — 20% at relevant distances. 

(ii) The overall good performance of the model does not mean that all results are 
good. For example, the model definitely 'overshoots' the repulsion in the isovector scalar 
mesonic channel | 19| , predicting a negative correlator in some window of distances. It 



violates positivity and is obviously wrong. 

(iii) Another sector, in which RIV does not perform correctly, is related to large-scale 
fluctuations of the topological charge. Obviously, a random distribution of instantons 
implies that in any volume V one expects to have excess of charge SQ ~ V 1 / 2 . In other 
words, the topological susceptibility x — < Q 2 > /V is non-zero. This feature is known to 
be wrong in the case of massless fermions, where \ should vanish for large volumes. Thus, 
one expects that the r\ correlator evaluated in the RIV 'overshoots' the right behaviour, 



which is indeed what was found 19 



The strong fluctuations of the topological charge can also be seen directly in the spec- 
trum of the Dirac operator and in the quark condensate. The generic argument for the 
eigenvalue distribution goes as follows. If the limit N c — > oo is taken before the thermo- 
dynamic limit, the matrix elements behave as independent random variables resulting in 
a semi-circular eigenvalue distribution fl5| , p6| . However, if the thermodynamic limit is 



taken at finite N c , as it should, the number of matrix elements is much larger than the 
number of independent collective variables, leading to strong correlations between them. 



In such case one expects the resulting eigenvalue distribution to be gaussian [13, On 
the other hand, an excess SQ of topological charge may create extra quasizero modes, but 
the naive argument^ giving rise to a divergent spectral level density p(A = 0) disagrees 
with our numerical simulations of the RIV. This is hardly surprising since the eigenvalues 
depend in a very nonlinear way on the matrix elements. 

Numerical results for an ensemble of 2560 configurations of 128 instantons (dotted line) 



5 If one ignores the interaction with the zero modes of the 'Q — background', one finds for a topological 
charge excess of SQ ~ \fV and corresponding overlap overlap matrix elements of T ~ p 2 /V 3 ^ 4 a level 
spacing of ~ leading to a divergent spectral density. 



and 160 configurations of 512 instantons (full line) are shown in Fig. 1. In both cases the 
density of instantons N/V = 1. The two spectra nearly coincide and are indistinguishable 
for small A, even if the bin size is decreased by an order of magnitude. This suggests that 
the thermodynamic limit will not alter our results. The dashed line is the corresponding 
gaussian distribution with variance given by a 2 = 2 < Tr(TT*)/N >, where < •■■ > 
denotes ensemble averaging. At small A our results suggest a strong departure from the 
gaussian distribution, which is related to the strong fluctuations of the topological charge. 
Although the spectrum itself does not become infinite, it seems to show an infinite slope 
p'(A) at A = 0. In terms of the normalized cumulants k^/k, 2 , and with values of 

0.180 and —0.910, respectively, the deviations from the gaussian distribution are much 
less visible. We also have checked that a semi-circular distribution is obtained when the 
overlap matrix elements are distributed as independent gaussian random variables with 
zero mean. 

The quark condensate, evaluated from the spectral density at A = 0, is enhanced by 
a significant factor, and fluctuations of the quark condensate, invoking the fluctuation- 
dissipation theorem, diverge in the chiral limit. This phenomenon should not take place 
if the inter-instanton interaction due to massless quarks is present. Indeed, simulations of 
an ensemble of interacting instantons do not show this enhancement at A = 0. A similar 
artificial peak at A = should be observed in quenched lattice calculations, but it should 
be absent in simulations with dynamical fermions. However, existing lattices are far too 
small to contain hundreds of instantons and to reveal this phenomenon clearly. 

In practice, the effect of this enhancement is small due to the finite quark masses 
which smear out the small eigenvalues. Considering the peak at A = as an artifact 
of the quenched approximation we will compensate for this by taking the quark masses 
somewhat larger than their phenomenological values, see details belowQ. 

3. Quark propagators: generalities 

In this section we discuss the quark propagator in a multi-instanton gauge field. Before 

6 In the case of dynamical quarks we have another reason for not taking the quark masses too small: 
at finite volume the smallest eigenvalues are ~ 1/V leading to a dip in the spectrum at A = 0. 
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going into details, let us clarify the following important point. Strictly speaking, the 
propagator is not a physical quantity because it is gauge-dependent. As discussed in 



detail in e.g. [fL7|, this problem may be circumvented d la Schwinger by complementing 
its definition with the path ordered exponential P exp[(ig/2) J% A a a^dx^)}. Physically, 
this corresponds to the addition of a static anti-quark: thus one naturally proceeds to the 
discussion of heavy- light mesons (see below). 

Although propagators are not physical quantities, the specific feature of the 'instan- 
ton vacuum' (which contains neither perturbative nor confining vacuum fields!) is that 
this extra Schwinger factor is actually very close to one: static quarks essentially ignore 
instantons. This statement was demonstrated explicitly in [|I7J], and a related analytical 
argument was given in [27]. In this sense, one may say that our data on propagators are 



'practically gauge invariant', and therefore are not purely academic quantities!]. 
The total spinor propagator 

S{x,y) =<x\ . - 1 . \y >, (3.1) 



will be approximated by the sum 



< x\i >< i\ . 1 \J>< J\y > +E < x \p >< p\ . n \ . \Q >< Q\y >, 

(3.2) 

where / and J are zero modes, or the 'bound states' in terms of the analogy put forward in 
chapter 2, and P and Q are nonzero modes which can be thought of as 'scattering states'. 
The latter can only be numerated by a continuum of eigenvalues, e.g. the momenta 
at large distances which makes the complete diagonalization of the Dirac operator a 
very complicated task. The simplifying assumption we have made is that the cross term 
between both type of modes is small and can be neglected. The first term follows from 
the Dirac operator in the space of zero modes (see eq. (2.5)) and is given by 

S ZM {x,y)=h(x)(-^—] <t>*j{y). (3.3) 

The second term is written as the sum of two contributions 

<x \—( F7vT~i ^il ■ D\y >' +im < x\— — — -\y >', (3.4) 



7 Readers not satisfied by this argument should note that the main objects of these papers, the corre- 
lation functions, are manifestly gauge invariant by construction. 
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where the ' indicates that the zero mode contribution should be excluded in the calcu- 
lation of the expectation value. In the second term we will neglect the interaction of the 
spin-induced gluomagnetic moment with the field. This amounts to the replacement of 
(7 • D) 2 — > D 2 , and is exact for some correlators [28]. What remains is the propagator of a 



scalar quark which will be discussed in section 4. The first term is the most complicated 
one and we postpone its discussion until section 5. 

The propagator in the approximation (3.1) still contains all possible gamma matrix 
structures. In other words, if we write 

S = $> r * (3.5) 
Ti = 1, 75, 7/» «757m> ply (p^v), ( 3 - 6 ) 

all coefficients Oj, which are S{7(3)-matrices, are nonzero. However, only two coefficients 
have a nonzero average value: a\ and a 7o (we assume that the propagation takes place in 
time direction). On the other hand, the average of the combination Tra^cif 7^ for all 
structures. In fact, linear combinations of these quantities are nothing else but mesonic 
correlation functions which will be discussed in [19| . 



The ensemble average of the propagator will be evaluated numerically via a Monte- 
Carlo simulation of 256 instantons in a box of 3.36 3 x 6.72/m 4 . The distribution of the 
position of the instantons will be taken uniform, whereas the size is kept fixed at 0.35/m. 
The orientations are sampled from the invariant group measure. The average propagators 
< S(x + r, x) > are calculated by averaging over an ensemble of 50 configurations and 
over 100 randomly chosen initial points x for each configuration and for each value of the 
separation r from the initial point. To u and d quark masses are taken equal to 10 MeV, 
which is still larger than the physical value. As discussed in section 2, the effect of the 
anomalous large number of small eigenvalues is suppressed this way (see Fig. 1). The 
strange mass is taken equal to the 'common sense' value of 140 MeV, 



4. Scalar quarks: propagators and mesons 

In this section we study the propagator of a massless scalar quark (or squark, for 
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brevity). Its analytical expression for one instanton is knownF] ||29|| : 



D(x,y) 



4tt 2 (x - y) 2 j\ + pi/( x - z )< 



[x — z) a 



(y 



(x — z) 2 (y — zy 



y/1 + p7( y - z y 



where z denotes the position of the center of the instanton. The propagator in the field 
of an anti- instanton is obtained by interchanging <j + and a~ . If the instanton is rotated 
in color space by the matrix R, the propagator should be rotated by the matrix R ab a b . 

The first factor Dq = l/47r 2 (x — y) 2 is the free scalar propagator. When all distances 
involved exceed the instanton size p, the remaining factors approach 1, leaving only the 
free propagator. Slightly less trivial is the short- distance limit (x — y) 2 — > 0. Expanding 
in powers of (x — y) 2 we obtain 

D , ^ = 1 + if^r}%,{x - z)„{y - z) v a a p 2 + _ _ _ 

' Att 2 (x — y) 2 4n 2 (x — y) 2 (x — z) 2 (y — z) 2 8tt 2 (x — z) 2 (y — z) 2 

and as the numerator of the second last term is proportional to (x — y), the free propagator 
Do dominates at small distances as well. 

The very complicated problem of propagation in a multi-instanton environment is 
significantly simplified at not too large distances, where only one or a few close instantons 
contribute. Analyzing the expression given above, one realizes that significant corrections 
only appear if either x or y or both are well within an instanton, i.e. at a distance no 
more than « p from its center. In |J only the effect of the 'closest' instanton was taken 
into account. 

One might think that one could improve on this by treating the effect of the other 
instantons as small perturbations and simply adding their contributions, 

D sc = D + J2(Dj~D ). (4.3) 



As discussed in [PH |, this correction is the first term of a 'rescattering' series. However, 



there is a serious problem: the integrated contribution of distant instantons diverges 
because of the term / d 4 z/(x — z) 2 {y — z) 2 . 



8 Here and below the instanton field is assumed to be in the singular gauge, A a ^ — 2r[^ v x v p 2 /[x 2 (x 2 +p 2 ], 
where rft v is the 't Hooft symbol. 
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In order to say more about of the origin of this term let us look at the expansion (4.2) 
in greater detail. When we allow for the possibility of a nonzero current squark mass the 
expansion of l/(—D 2 + m 2 ) in powers of the gauge field results in 

_ n2 \ — ^(x,y) = D m (x,y) + f d 4 uD m (x,u)ia a A a ^u - z)d fl D m (u,y) + ■ ■ ■ (4.4) 

This expansion which also was considered in |30j provides us with the corrections of 0(p 2 ). 
The interpretation of the correction in this formula is that the quark propagates to some 
point u in space time, then exchanges 2 or more gluons with the center of the instanton 
(obtained from the expansion of the instanton profile in inverse powers of (u — z) 2 ) and 
finally proceeds to its end point. Therefore the factors l/(x — z) 2 and l/(x — y) 2 , obtained 
after integration over u do not have an unambiguous interpretation as the propagator of 
a gluon or of a squark. Both gluons and squarks acquire a mass due to non-perturbative 
effects, and therefore do not propagate far. Writing the propagator as 

Dfay) = flb(»,») + zy SD ^_ z)2 , (4.5) 

suggests a cure to our 'infrared problem': one should substitute the massless propagators 
by the massive propagators Dm{x, z) and Dm(v, z), which makes the integral convergent^]. 
Since the divergence is logarithmic, the exact form of the cutoff is not very important. 
The value of the mass is certainly larger than the current quark masses. The natural 
scale for this cut-off mass M ~ Aq C d- In the case of the strange quark one might argue 
to add its mass to the cutoff mass, but since we are already considering a relative small 
correction, we will refrain from this kind of fine tuning. 

Results for TrD(x,0) as a function of the distance are shown in Fig. 2a. For con- 
venience, TtD(x, 0) is normalized to the free massless scalar propagator D = l/47r 2 x 2 , 
so the deviation from 1 is caused by the deflection of the squark in the instanton gauge 
field. The solid line corresponds to a massive scalar propagator, with a fitted mass value 
of M squark = 140 MeV. Note, that this value is roughly half the 'constituent quark mass', 
expected for spinor quarks from phenomenological considerations. 

In Fig. 2b we show the simplest mesonic correlation function, obtained by averaging 



9 This solution to our infrared problem is in fact similar to the well known diagram resummation in 
plasma physics which takes care of Debye screening. 
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the square of the squark propagator 

K(x) =< Tr(D(x, 0)D(0, x)) > . (4.6) 

As usual, the data shown are normalized with respect to the perturbative correlator, and 
should be close to 1 at small x. The solid line corresponds to the square of the massive 
propagator, fitted to data in Fig. 2a. If two quarks would propagate independently, 
the data would be in agreement with it. However, all measured points are above this 
curve, which means that scattering on instantons not only generates an effective squark 
mass, but also an attractive interaction. The fitted value of the meson mass is ~ 150MeV, 
which is less than 2M squark rs 280MeV. Of course, it is natural that due to such attractive 
interaction scalar quarks form mesons. 

Note, that the origin of this attraction at small distances can be traced to the second 
term in the r.h.s. of (4.2): its average (e.g. over orientations) is zero, but not the average 
of its square. It is the simplest example of 'hidden' components of the quark propagators 



to be discussed in the next paper of this series |L9 . 

Concluding this section, we have shown that the instanton vacuum can produce an 
effective squark mass and an interaction between them. This implies that similar effects 
can be expected for spinor quarks as well. 



5. Effect of non-zero modes for spinor quarks 

The main special feature of the spinor propagator is that it also receives a contribution 
from the zero modes. However, before discussing the complete propagator we first study 
the part of non-zero mode contribution that is given by the first term in eq. (3.4). It is 
known analytically for a massless quark in the field of a single instanton [^9] : 

ry I \ 1 1 ( n I \ /-, p 2 <J~ ■ X <7 + ■ y , 

Si{x,y) = . , \S {x,y){l + — ) 

j l + p 2 /x 2 Jl + p 2 /y 2 V X V 



p 2 a ■xa + -~/a -Aa + -y + a -xa + -Aa ■ 7 cr + ■ y \ 

p 2 + x 2 75 + p 2 + y 2 75 7 ' 

(5.1) 



where we have introduced the projectors 75 = (1 ± 7s)/2. The massless free quark 
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propagator is denoted by So(x,y), and the free massless scalar propagator, D (x,y) is 
given by the first term in the r.h.s. of eq. (4.6). The expression for anti-instantons 
is obtained by interchanging a + and a~ . The SU(3)— extension of the propagator is 
obtained by embedding this propagator in the upper 2x2 block of a 3 x 3 matrix, and 
substituting the free propagator in the remaining diagonal matrix element. For arbitrary 
orientation, U, of the instanton Sj(x,y) has to be rotated covariantly. 

We generalize Sj(x,y) to many instantons, in the same approximate way as in the 
previous section for squarks, namely by summing all deviations from the free propagator: 



Although in this case the contribution of distant instantons is not divergent, we still think 
it should be excluded on the basis that gluons cannot propagate that far. If the NZM- 
induced effects are treated as a correction, it is logical to also include here an effective 
mass to cut off the contribution of distant instantons. 

The last necessary step is to include the effects of a non-zero bare quark mass: this 
is especially needed for the discussion of correlators involving strange quarks. Generally 
speaking, there are at least three different scales to which the strange quark mass, m s , 
should be compared: 

(i) The eigenvalues of the Dirac operator in the ZMZ which are comparable to m s . There- 
fore, one should be careful at this point and not treat m s as a small parameter^]: the 
mass in the strange component of the denominator of the ZMZ term in the propagator 
will be taken equal to m s . 



(ii) The field strength of the gauge field inside instantons. Its value is very large: JG 2 ^ ~ 

8\^3/gp 2 ~ 2GeV 2 ^> m 2 ~ .02GeV 2 , and therefore the motion of strange quarks inside 
instantons is essentially the same as that of massless ones. 

(iii) The total (Euclidean) propagation time. In this work this time is comparable to the 
strange quark mass (r ~ 1 — 2fm ~ w-^ 1 ), which makes the strange quark propagators 
substantially different from non-strange ones. However, in a reasonably dilute instanton 
vacuum most of this motion takes place in relatively empty space, and one may therefore 
describe these corrections using the free massive propagator. 

10 Note that at this point we deviate from what is done in chiral perturbation theory. 




(5.2) 
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In such a complex situation, the problem is to work out a useful interpolating for- 
mula, which combines several limiting cases for which an analytical solution is known and 
presumably provides a reasonable description in the general case. After trying several ex- 
pressions, we arrived at the following expression in which factorization of the mass-related 
damping of the propagator from corrections of scattering on instantons is assumed: 

5 7 (m, x, y) = S m (x, y) + D M (x, 0)6 S(x, y)D M (0, y). (5.3) 

The correction SS is given by 

D o {x,0)D o {p,y) 

and S m (x, y) is the massive free fermion propagator. At small distances (inside instantons) 
the the propagation of quarks approaches free propagation, while at large \x—y\ the quarks 
move as free massive fermions in between successive scatterings on instantons while the 
correction 6S is exponentially suppressed. We will choose M = Aqcd which we believe 
to be a natural choice for this cut-off mass. 

Finally we proceed to the expression for the total propagator S(x, y) by adding the 
contributions discussed is this section and in section 3 and 4. As a result we find 

S(x, y) = 5 7 (m, x, y) + (j>i{x) ( , ) <f>j(y) + imD sc (x, y), (5.5) 

where D sc is defined in eqs. (4.3) and (4.5) modified as discussed below eq. (4.5). 

In Fig. 3 we show the chirality flip Tr S (a) and the chirality non-flip Tr 70S' (b) com- 
ponents of this propagator. The squares represent results for the full propagator, whereas 
the crosses are for the propagator S(x, y) = S m (x, y) + S ZM (x, y), which includes only the 
modification of the propagator due to zero modes. The first trace Tr S is normalized with 
respect to the short distance limit of the massive free quark propagator Tr So (x,m)/m 
such that its value at x = is equal to the mass m. The second trace is normalized with 
respect to the free propagator. One observes that their effect is significant, especially at 
medium distances x ~ 1/2 fm. 

The physical question we are going to discuss is whether or not these results can be 
interpreted as the appearance of some universal effective mass m e ^/Q For this reason we 

11 The existence of a momentum dependent effective mass has also been obtained analytically in the 
large N c - limit |]|. 
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have plotted the behaviour of the massive fermion propagator for m = 200, 300, 400 MeV. 
At small x, the calculated points start from zero, consistent with bare quark propagation 
with a small input mass. However, for x > 0.5/m they even somewhat 'overshoot' the 
curves, but qualitatively imitate their x— dependence. Since none of the curves is par- 
ticularly close to the data trend it is not possible to select a 'preferable' mass from the 
chirality flip amplitude. 

The chirality non-flip component of the propagator, shown in Fig. 3b, tends to 1 at 
small x (which corresponds to massless propagation), but at larger x we do find qualitative 
agreement with the massive free propagators. The square points follow the pattern given 
by the solid line. This means that the net effect of the 'instanton vacuum' on quark 
propagation can indeed be represented approximately by the appearance of an effective 
mass of the order of 300 MeV. This is about twice larger than for the scalar quark. 
Indeed, when only the contribution of the zero modes is included, as shown by the crosses 
in Fig. 3b, the 'effective mass' is about twice smaller. This shows that roughly half the 
constituent mass is due to the zero modes and the other half due to the non-zero modes. 

Note that the curves in Fig. 3 imply that the amplitude of 'bare-to-dressed' quark 
transition amplitude (the Z-factors introduced in ||16|| ) is exactly 1, so this feature of TOO 
percent dressing' is approximately reproduced by our propagators as well. 

In conclusion, there is some qualitative agreement with the 'constituent quark' model. 
However, one should not take it literally. First of all, we remind that agreement is not 
really good for the chirality-flip amplitude shown in Fig. 3a. And, last but not least, if two 
(or more) quarks propagate together, they show a very strong interaction. Therefore, as 
we will show in future publications, none of the mesonic or baryonic correlators measured 
are actually reproduced by the naive 'constituent quark' model. 



6. Heavy-light mesons 

Heavy-light mesonic correlators were first discussed by one of us in ||32|| , where the 
concept of 'heavy quark symmetry' between any channels with sufficiently heavy quark 
masses was introduced. Indeed, in the heavy quark limit Mq ^> Aqcd the heavy quark 
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behaves as a static center (the mesons' 'nucleus'), and the corresponding correlation func- 
tions essentially reduce to the propagator of the light quark. As the direction of the spin 
of heavy quark is irrelevant in this static limit, a specific degeneracy of correlation func- 
tions follows. For example, the spatial component of the vector correlator coincides with 
pseudoscalar correlator, and the spatial component of the axial correlator coincides with 
scalar correlator. 

More specifically, we discuss correlators of currents containing one heavy Q = (c, b, t) 
and one light quark q defined by the expectation value 

K r (x) =< QTqqTQ >, (6.1) 

where V is one of the possible gamma matrix structures. In this section we consider 1 
(S), 75 (P), 7 (V) and 775 (A). The heavy quark propagator that enters (6.1) is given 
by the nonrelativistic quark propagator times a path ordered exponential. In |l7j it was 
shown that the correction due to the path ordered exponential is small. This allows us 
to translate the propagator immediately into heavy-light correlation functions. Using 
that only the large components of the heavy quark spinor are relevant one obtains the 
correlator 

K T (x) = zTr(ff(-x)r 1+7 ° x/|x| r) (6.2) 

times a non-relativistic quark propagator. The separation x is chosen along the positive 
time axis. Only the parity of the current is essential, and for the parity P = ±1 channels 
with T = 1, 75 one obtains the following correlation functions 

K±(x) = tTr[(^^)S(-x)}. (6.3) 



As was noticed in [32|, the splitting of these correlators at small distances, normalized to 



the free quark correlator Kq, is simply proportional to the quark condensate 

2 3 

K ± (x)/K (x) = 1 T —ttH <4>ip >\-\ ■ (6.4) 

6 

In order to describe the correlator in the whole region, we use a standard parametrization 
for the spectral function, a resonance plus the perturbative continuum above a certain 
'threshold' energy E , in this case given by f32"f : 

lmK ± (E) = 6nn5(E - E ICS ) + 6(E - E )^, (6.5) 

Z7l 



where n = fqM q / 12 is the 3- dimensional density of the light quark at the center. In terms 
of the space-time correlator, it translates into the following expression: 

K±(x)/K (x) = 2Tx 2 nx 3 e~ E ^ x + (1 + E x + E 2 x 2 /2)e~ Eox . (6.6) 

Our results for the propagator, in the form K±(x)/Kq(x), are shown in Fig. 4, where 
the symbols V and A refer to the spatial components of the vector and axial currents, 
respectively. The dashed curves correspond to the three-parameter fit (6.6), with the 
values of the parameters given in Table 1. For comparison we also give results of previous 
works using various methods^- 

The resonance position can also be compared to the experimental values for B meson 
masses, provided that the b quark mass value is obtained from different sources (e.g. sum 
rules for upsilons). Assuming that m& pa 4800 MeV, as follows from such sources, one 
obtains for E ves = itlb — pa A75MeV, which is not too far from our fitted value of 615 
MeV. 

The splitting of opposite parity states is only known for charmed mesons (such states 
have not yet been discovered for B mesons): it is about 450 MeV (see discussion is |l|), to 
be compared to our estimated difference of 575 MeV . This value is very close to twice the 
constituent quark mass, a result which also was derived in using completely different 
methods. 

It is instructive to consider the small distance expansion of the correlation function 
(6.6), 

K(x)/K (x) = 1 + (2n 2 n - \e 3 )x 3 + {\e a - 2n 2 nE ms )x A + ■■■. (6.7) 

o 

Equating this to (6.4) we find the so called 'duality relation', 

±^<-W>= (2vr 2 n ± - \Ef), (6.8) 
6 6 

which should be approximately satisfied. From the r.h.s. we find the value for the quark 
condensate of — (215 MeV) 3 in case of the 0~ state and — (578 MeV) 3 in case of the + 
state. The discrepancy between these two numbers is a signature of an early deviation 
of our fit and, possibly, the true correlator from the predictions of the operator product 



12 Comparison of compiled lattice results with results from QCD sum rules, as well as the discussion 
of 1/rriQ corrections can be found e.g. in 
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expansion. This also follows from the right hand side of the dispersion relation which can 
be expanded in powers of x for x <C l/E ~ 0.2 fm only. 

7. Conclusions and discussion 

In this work we report an extensive numerical study of the simplest ensemble of instan- 
tons: the random instanton vacuum. This vacuum is characterized by two parameters: 
(i) the size p and (ii) the density N/V of instantons (plus anti-instantons) in the QCD 
vacuum, which were chosen to reproduce the size of the quark condensate and the gluon 
condensate. 

The topological charge density in this model fluctuates strongly on all length scales, 
which shows up in the form of an excess of small eigenvalues of the Dirac operator above 
the expected gaussian distribution. 

For simplicity, we first studied the propagation of scalar quarks (or squarks), in which 
case there is only a nonzero mode contribution to the propagator. It appears that the 
scalar mesonic correlation function can be described remarkably well by a 'constituent 
quark' model with a mass of about 140 MeV . 

The spinor quark propagator has two different components: the chirality non-flip 
and the chirality flip parts. The first one can again be described by an effective quark 
mass of m e ff ~ 300 MeV. The second component is more complicated and corresponds 
approximately to a free massive propagator at distances larger than about 1/2 fm. At 
distances of x ~ 0.5 — 1.0/m both the zero and non-zero fermion modes contribute a 
comparable amount to the constituent quark mass. 

Observable consequences of these statements are most clearly seen in the spectra of 
heavy-light mesons. We find a mass difference of 615 MeV between the B— meson and the 
bottom quark which agrees well with the phenomenological value. The splitting between 
the mass of opposite parity states was found to be 575 MeV, which is roughly equal to 
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twice the constituent quark mass. 
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channel J F 



E TCS [MeV] 



n 



f 2 M Q /12, [fm- 



E [MeV] 



ref. 



comment 



o-r 



615 

400±100 

480±80 

500±30 



1.2 

1±.5 

1±.3 

2. ±.5 

4.5±1 

2.5±1. 



1010 

900±100 

700 

1100±100 



this work 



m 

\M> HI 



random instantons 
QCD sum rules 
IIA 

sum rules 
lattice, 'static' 
lattice, 'dynamic' 



o + ,r 



1190 

1200±200 
1000± 200 



4.6 
6.3 
.7±.5 



1745 

1800±200 
1200±200 



this work 



random instantons 
QCD sum rules 
IIA 



Table 1. 



Fitted parameters of heavy-light mesons, compared to those derived in previous works. 
E res , n, E are the resonance mass (counted from the heavy quark mass), the density of 
the light quark at the center and the 'asymptotic freedom threshold' (see text). For 
comparison, we give some lattice results based on a static heavy quark and a 1/m.Q 
expansion, as well as results based on the extrapolation of simulations with dynamical 
quarks. 
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Figure Captions. 



Fig 1. The eigenvalue density n(A) of the Dirac operator in the space of zero modes for 
an ensemble of 512 and 128 instantons. The normalization Pf^ d\p(X) = 1 and the bin 
size used to obtain the histograms was 0.025. The variance of the gaussian curve is equal 
to the variance of n(A) for iV = 512. 

Fig. 2. The average propagator of a scalar quark (a) and its averaged square (b) (cor- 
responding to the scalar meson), normalized to the propagator of a massless quark. The 
solid lines in (a) and (b) correspond to a propagator with a fitted mass value of 140 MeV. 

Fig. 3. The chirality-flip (a) and non-flip (b) components of the quark propagator versus 
the distance x (in fm). The normalization, indicated in the figure, is discussed in the 
text. Crosses in (b) correspond to the effects of zero- modes only, while the squares give 
the complete result. Three lines, the short-dashed, solid and long-dashed correspond to 
a massive free propagator with a mass of 200,300 and 400 MeV, respectively. 

Fig. 4. The correlation functions for negative (P = pseudoscalar, V = vector (spatial)) 
and positive (S = scalar, A = axial (spatial)) parity heavy-light mesons. The data points 
are our results, while the dashed curves are fits, discussed in the text and in Table 1. 
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